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Abstract. We present a summary of the physics of gluon saturation and non-linear QCD evolution at small values 
of parton momentum fraction x in the proton and nucleus in the context of recent experimental results at HERA and 
RHIC. The rich physics potential of low-x studies at the LHC, especially in the forward region, is discussed and some 
benchmark measurements in pp, pA and AA collisions are introduced. 
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1 Introduction 

1 .1 Parton structure and evolution 

The partonic structure of the proton (nucleus) can be probed 
with high precision in deep inelastic scattering (DIS) electron- 
proton ep (electron-nucleus, eA) collisions. The inclusive DIS 
hadron cross section, d 2 o/dxdQ 2 , is a function of the virtuality 
Q 2 of the exchanged gauge boson (i.e. its "resolving power"), 
and the Bjorken-x fraction of the total nucleon momentum car- 
ried by the struck parton. The differential cross section for the 
neutral-current (y, Z exchange) process can be written in terms 
of the target structure functions as 



2%a 2 



dxdQ 2 xQ 4 



[Y+-F 2 TY_-xF 3 - 



-y 2 -F L ] 



(1) 



where Y± = 1 ± (1 — y) 2 is related to the collision inelasticity 
y, and the structure functions F 2 t,_l(x,Q 2 ) describe the density 
of quarks and gluons in the hadron: F 2 « e 2 x + #,-), xF-$ oc 
xT,j(qj — qi), Fi OL s xg {xqi,xqi and xg are the corresponding 
parton distribution functions, PDF). F 2 , the dominant contribu- 
tion to the cross section over most of phase space, is seen to 
rise strongly for decreasing Bjorken-x at HERA (Fig. 0. The 

growth in F2 is well described by F 2 (x, Q 2 ) x~^®~\ with X sa 
0. 1 - 0.3 logarithmically rising with Q 2 1 1 1. The F 2 scaling vio- 
lations evident at small x in Fig.^are indicative of the increas- 
ing gluon radiation from sea quarks. The xg(x, Q 2 ) distribution 
itself can be indirectly determined (Fig.|2ji from the F 2 slope: 



dF 2 ( Xl Q 2 ) _ 10a,(g 2 ) 
3ln(g 2 ) ~ 277t 



xg(x,Q 2 ). 



(2) 



Although the PDFs are non-perturbative objects obtained from 
fits to the DIS data, once measured at an input scale Qq > 
2 GeV 2 their value at any other Q 2 can be determined with 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evo- 
lution equations which govern the probability of parton branch- 
ings (gluon splitting, q,g— strahlung) in QCD [3 1. 



HERA F 




Q 2 (GeV 2 ) 

Fig. 1. F 2 (x, Q 2 ) measured in proton DIS at HERA ( s/s = 320 GeV) 
and fixed-target ( *Js ~ 10-30 GeV) experiments 

The DGLAP parton evolution, however, only takes into ac- 
count the (2 2 -dependence of the PDFs, effectively summing 
leading powers of [a s ln((2 2 )] ("leading twist") generated by 
parton cascades in a region of phase space where the gluons 
have strongly-ordered transverse momenta towards the hard 
subcollision Q 2 3> k 2 T ^> ■ ■ ■ ^> k 2 T . Such a resummation is 
appropriate when ln(2 2 ) is much larger than ln(l fx). For de- 
creasing x, the probability of emitting an extra gluon increases 
as oc a s ln(l /x). In this regime, the evolution of parton densi- 
ties proceeds over a large rapidity region, Ay ~ ln(l/x), and 
the finite transverse momenta of the partons become increas- 
ingly important. Here the full kj phase space of the gluons (in- 
cluding scattering of off-shell partons) has to be taken into ac- 
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Fig. 2. Gluon distributions extracted at HERA (HI and ZEUS) as a 
function of x in three bins of Q 2 |5) 

count and not just the strongly-ordered DGLAP part. Thus, the 
appropriate description of the parton distributions is in terms 

of kj-unintegrated PDFs, xg(x 1 Q 2 ) = f@ dkjG(x,kj) with 
G(x,kj,) ~ h(kj)x~^, described by the Balitski-Fadin-Kuraev- 
Lipatov (BFKL) equation |4| which governs parton evolution 
in x at fixed Q 2 . Hints of extra BFKL radiation have been 
recently found at HERA in the enhanced production of for- 
ward jets compared to DGLAP expectations 15 161 . At large Q 2 , 
a description resumming over both a s ln(2 2 ) and a$ln(l/x) 
is given by the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) 
evolution equation j7| . 

1.2 Parton saturation and non-linear evolution at Iow a 

As shown in Figure the gluon density rises very fast for 
decreasing x. Eventually, at some small enough value of x 
(a. s ln(l/x) ^> 1) the number of gluons is so large that non- 
linear (gg fusion) effects become important, taming the growth 
of the parton densities. In such a high-gluon density regime 
three things are expected to occur: (i) the standard DGLAP and 
BFKL linear equations should no longer be applicable since 
they only account for single parton branchings (1 — >2 pro- 
cesses) but not for non-linear (2 — >1) gluon recombinations; 
(ii) pQCD (collinear and kj) factorization should break due to 
its (now invalid) assumption of incoherent parton scattering; 
and, as a result, (iii) standard pQCD calculations lead to a vio- 
lation ofunitarity even for Q 2 ^> Ag CD . Figure[3]schematicaliy 
depicts the different parton evolution regimes as a function of 
y — \n(l/x) and Q 2 . For small enough x values and for vir- 
tualities below an energy-dependent "saturation momentum", 
Q s , intrinsic to the size of the hadron, one expects to enter the 
regime of saturated PDFs. Since xg(x,Q 2 ) can be interpreted 
as the number of gluons with transverse area r 2 ~ 1 /Q 2 in the 
hadron wavefunction, an increase of Q 2 effectively diminishes 
the 'size' of each parton, partially compensating for the growth 
in their number (i.e. the higher Q 2 is, the smaller the x at which 




Fig. 3. QCD "phase diagram" in the 1 /x, Q 2 plane (each dot repre- 
sents a parton with transverse area ~ 1 /Q 2 and fraction x of the hadron 
momentum). The different evolution regimes (DGLAP, BFKL, satura- 
tion) as well as the "saturation scale" and "geometric scaling" curves 
between the dense and dilute domains are indicated. Adapted from 1 8 1 



saturation sets in). Saturation effects are, thus, expected to oc- 
cur when the size occupied by the partons becomes similar to 
the size of the hadron, %R 2 . In the case of nuclear targets with 
A nucleons (i.e. with gluon density xG = A -xg), this condition 
provides a definition for the saturation scale 191 101 : 

Q 2 (x) ~ a s ^xG(x,Q 2 ) ~ A 1 ^,-* ~A x /3( v^A^V*, 

(3) 

with A, f» 0.25 1 1 i jj. Eq. tell us that Q s grows with the num- 
ber of nucleons in the target and with the energy of the colli- 
sion, -s/s, or equivalently, the rapidity of the gluon y = ln( 1 jx) . 
The nucleon number dependence implies that, at equivalent en- 
ergies, saturation effects will be enhanced by factors as large as 
A 1 / 3 w 6 in heavy nuclear targets (A = 208 for Pb) compared 
to protons. In the last fifteen years, an effective field theory 
of QCD in the high-energy (high density, small x) limit has 
been developed - the Colour Glass Condensate (CGC) - 
which describes the hadrons in terms of classical fields (satu- 
rated gluon wavefunctions) below the saturation scale Q s . The 
saturation momentum Q s introduces a (semi-)hard scale, Q s ^> 
Agco, which not only acts as an infrared cut-off to unitarize 
the cross sections but allows weak-coupling perturbative cal- 
culations (OL s (Q s ) <Cl) in a strong colour field background. 
In the CGC framework, hadronic and nuclear collisions are 
seen as collisions of classical wavefunctions which "resum" all 
gluon recombinations and multiple scatterings. The quantum 
evolution in the CGC approach is given by the JIMWLK 1 13 1 
non-linear equations (or by their mean-field limit for N c — > °°, 
the Balitsky-Kovchegov equation 1141 ) which reduce to the 
standard BFKL kernel at higher x values. 



2 Parton saturation: Experimental studies 

The main source of information on the PDFs is obtained from 
hard processes as they involve outgoing particles directly cou- 
pled to the partonic scattering vertices. Figure |4] summarizes 
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Fig. 4. Examples of experimental measurements at various facilities 
providing information on the gluon PDF in the range x ~ 10~ 5 — 0.8 



the variety of measurements at different experimental facilities 
which are sensitive to the gluon density and their approximate 
x coverage. xG enters directly at LO in hadron-hadron colli- 
sions with (i) prompt photons, (ii) jets, and (iii) heavy-quarks 
in the final state, as well as in the (difficult) DIS measurement 1 
of (iv) the longitudinal structure function Fl- In addition, (iv) 
heavy vector mesons T) from diffractive photoproduction 
processes 2 are a valuable probe of the gluon density since their 
cross sections are proportional to the square of xG | 1 61171 : 



da 



ypA^V pA 



dt 



3aMl 



■167t 3 [xG(x,Q 2 )] 2 , 



with Q 2 =M v /4 and x = M v /Wy p A . 



(4) 



(5) 



practice the (xi) values in parton-parton scatterings are at least 
10 larger than x™" 1181 ). it provides the right estimate of the 
typical X2 = (pr/ ^ values reached in non-linear 2 — ► 1 

processes (in which the momentum is balanced by the gluon 
"medium") as described in parton saturation models 1 1 91201 . 

Figure[5]shows the kinematical map in (x, Q 2 ) of the exist- 
ing DIS, DY, direct y and jet data used in the PDF fits. Results 
from HERA and the Tevatron cover a substantial range of the 
proton structure (10~ 4 < x < 0.8, 1 < Q 2 < 10 5 GeV 2 ) but the 
available measurements are much rarer in the case of nuclear 
targets (basically limited to fixed-target studies, 10~ 2 < x < 0.8 
and 1 < Q 2 < 10 2 GeV 2 ). As a matter of fact, the nuclear parton 
distributions are basically unknown at low x (x < 0.01) where 
the only available measurements are fixed-target data in the 
non-perturbative range (Q 2 < 1 GeV 2 ) dominated by Regge 
dynamics rather than quark/gluon degrees of freedom. An ex- 
ample of the current lack of knowledge of the nuclear densi- 
ties at low x is presented in Fig. |6] where different available 
parametrizations of the ratio of Pb to proton gluon distribu- 
tions, consistent with the available nDIS data at higher x, show 
differences as large as a factor of three at x ~ 10~ 4 12 3 1241 . 




The main source of information on the quark densities is ob- 
tained from measurements of (i) the structure functions in 
lepton-hadron scattering, and (ii) lepton pair (Drell-Yan, DY) 
production in hadron-hadron collisions. In hadronic collisions, 
one commonly measures the perturbative probes at central ra- 
pidities (y = 0) where x=xj = QJ ^/s, and Q ~ pr,M is the 
characteristic scale of the hard scattering. However, one can 
probe smaller xi values in the target by measuring the corre- 
sponding cross sections in the forward direction. Indeed, for 
a 2 — > 2 parton scattering the minimum momentum fraction 
probed in a process with a particle of momentum pj produced 
at pseudo-rapidity r| is [ 18 1 
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where XT—2pr/y/s, 



(6) 



i.e. x% m decreases by a factor of ~10 every 2 units of rapidity. 
Though Eq. is a lower limit at the end of phase-space (in 

xG can also be (in)directly extracted from Fi through the deriva- 
tive in Eq. J5J as well as from the f^ harm data [ l5l . 

2 Diffractive yp (yA) processes are characterized by a quasi-elastic 
interaction - mediated by a Pomeron or two gluons in a colour singlet 
state - in which the p (A) remains intact (or in a low excited state) and 
separated by a rapidity gap from the rest of final-state particles. 



Fig. 5. Available measurements in the (x, Q 2 ) plane used for the deter- 
mination of the proton |21 1 (top) and nuclear |22| (bottom) PDFs 



2.1 HERA results 

Though the large majority of ep DIS data collected during 
the HERA-I phase are successfully reproduced by standard 
DGLAP predictions, more detailed and advanced experimen- 
tal and theoretical results in the recent years have pointed to 
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Fig. 6. Ratios of the Pb over proton gluon PDFs versus x from different 
models at Q 2 = 5 GeV 2 . Figure taken from 1231 

interesting hints of non-linear QCD effects in the data. Ar- 
guably, the strongest indication of such effects is given by the 
so-called "geometric scaling" property observed in inclusive 
Odis for* < 0.01 1|25 1 as well as in various diffractive cross sec- 
tions 1261271 . For inclusive DIS events, this feature manifests 
itself in a total cross section at small x (x < 0.01) which is only 
a function of x = Q 2 /Q 2 S (x), instead of being a function of x and 
Q 2 /Q 2 separately (Fig.0. The saturation momentum follows 
Qs{x) — Qaix/xf))^ with parameters A. ~ 0.3, Qo = 1 GeV, and 
xq ~ 3T0~ 4 . Interestingly, the scaling is valid up to very large 
values of x, well above the saturation scale, in an "extended 
scaling" region (see Fig.|3 where Q 2 < Q 2 < Q 4 JA 2 QCD [28 
8|. The saturation formulation is suitable to describe not only 
inclusive DIS, but also inclusive diffraction y*p — > X p. The 
very similar energy dependence of the inclusive diffractive and 
total cross sections in y*p collisions at a given Q 2 is easily ex- 
plained in the Golec-Biernat Wiisthoff model 1 25 1 but not in 
standard collinear factorization. Furthermore, geometric scal- 
ing has been also found in different diffractive DIS cross sec- 
tions (inclusive, vector mesons, deeply-virtual Compton scat- 
tering DVCS) I26I27I . All the observed scalings are suggestive 
manifestations of the QCD saturation regime. Unfortunately, 
the values of Q 2 ~ 0.5 - 1.0 GeV at HERA lie in the transi- 
tion region between the soft and hard sectors and, therefore, 
non-perturbative effects obscure the obtention of clearcut ex- 
perimental signatures. 



2.2 RHIC results 

The expectation of enhanced parton saturation effects in the nu- 
clear wave functions accelerated at ultra-relativistic energies, 
Eq. 0, has been one of the primary physics motivations for the 
heavy-ion programme at RHIC 3 1 12 1. Furthermore, the proper- 
ties of the high-density matter produced in the final-state of 
AA interactions cannot be properly interpreted without having 
isolated first the influence of initial state modifications of the 

3 At y =0, the saturation scale for a Au nucleus at RHIC (Q 2 ~ 2 
GeV) is larger than that of protons probed at HERA (Q 2 S ~ 0.5 GeV). 
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Fig. 7. Geometric scaling in the DIS y*p cross sections plotted versus 
X = Q 2 /Q 2 in the range x < 0.01, 0.045 < Q 2 < 450 GeV 2 ED 

nuclear PDFs. After five years of operation, two main experi- 
mental results at RHIC have been found consistent with CGC 
predictions: (i) the modest hadron multiplicities measured in 
AuAu reactions, and (ii) the suppressed hadron yield observed 
at forward rapidities in dAu collisions. [In addition, a recent 
analysis of the existing nuclear DIS F2 data also confirms the 
existence of "geometrical scaling" for* <0.017 |29|.] 

Vs (GeV/A) 




Fig. 8. Dependences on cm. energy and centrality (given in terms 
of the number of nucleons participating in the collision, A^ part ) of 
rfA , c /,/dr|| T1= o (normalized by A'part): PHOBOS AuAu data 1 34| vs the 
predictions of the saturation approach 1 29 

The bulk hadron multiplicities measured at mid-rapidity 
in central AuAu at i/s^ = 200 GeV are dN c h / d^^o w 700, 
comparatively lower than the dN c f, /c/r| |n=o ~ 1000 expecta- 
tions of "minijet" dominated scenarios |30|, soft Regge mod- 
els 1 3 1 1 (without accounting for strong shadowing effects 1 32 1), 
or extrapolations from an incoherent sum of proton-proton col- 
lisions 1 33 1 . On the other hand, CGC approaches 1 1 1 29 1 which 
effectively take into account a reduced number of scattering 
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centers in the nuclear PDFs, f a /A(x, Q 2 ) < A ■ f a /N{ x >Q 2 ) re- 
produce well not only the measured hadron multiplicities but - 
based on the general expression Q - also the centrality and 
cm. energy dependences of the bulk AA hadron production 
(Fig.®. 



10" 



d+Au -» h +X @ v s NN = 200 GeV (ti = 3.2) 

>V BRAHMS data 

— pQCD+shadow [Guzey et al.] 

— pQCD+shadow [Accardi] 

CGC [Tuchin et al.] 

CGC [Jalilian-Marian] 




4 5 
P T (GeV/c) 

Fig. 9. Nuclear modification factor RcIAu{pt) for negative hadrons at 
T| =3.2 in dAu at = 200 GeV: BRAHMS dat a l35l c ompared to 
leading-twist shadowing pQCD 1 18 19 1 and CGC 1 36 37 1 predictions 

The second manifestation of saturation-like effects in the 
RHIC data is the BRAHMS observation (35) of suppressed 
yields of moderately high-pr hadrons (pj « 2 — 4 GeV/c) in 
dAu relative to pp collisions at r| ps 3.2 (Fig. [9). Hadron pro- 
duction at such small angles is sensitive to partons in the Au 
nucleus with X2 ~ <??(10~ 3 ). The observed nuclear modification 
factor, RdAu ~ 0.8, cannot be reproduced by pQCD calculations 
that include standard leading-twist shadowing of the nuclear 
PDFs 1181191 but can be described by CGC approaches 1361371 
that parametrize the Au nucleus as a saturated gluon wavefunc- 
tion. As in the HERA case, it is worth noting however that at 
RHIC energies the saturation scale is in the transition between 
the soft and hard regimes (Q 2 2 GeV 2 ) and the results consis- 
tent with the CGC predictions are in a kinematic range with rel- 
atively low momentum scales {(pr) ~ 0.5 GeV for the hadron 
multiplicities and (pj) ~ 2.5 GeV for forward inclusive hadron 
production) where non-perturbative effects can blur a simple 
interpretation based on partonic degrees of freedom alone. 



3 Low-jc QCD at the LHC 

The Large Hadron Collider (LHC) at CERN will provide pp, 
pA and AA collisions at ^/~s~^ = 14, 8.8 and 5.5 TeV respec- 
tively with luminosities Jzf ~ 10 34 , 10 29 and 5-10 26 cirr 2 s -1 . 
Such large cm. energies and luminosities will allow detailed 
QCD studies at unprecedented low x values thanks to the copi- 
ous production of hard probes (jets, quarkonia, heavy-quarks, 
prompt y, Drell-Yan pairs, etc.). The advance in the study of 
low-x QCD phenomena will be specially substantial for nuclear 
systems since the saturation momentum, Eq. (|3}, Q 2 « 5 - 10 
GeV 2 , will be in the perturbative range 1111 . and the relevant x 
values, Eq. l|6}, will be 30-70 times lower than AA and pA re- 
actions at RHIC: x w 10 3 (10 5 ) at central (forward) rapidities 
for processes with Q 2 -10 GeV 2 (Fig.fTol. 
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Fig. 10. Kinematical (x, Q 2 ) range probed at various rapidities y and 
cm. energies in ,y5W = 8.8 TeV pA collisions at the LHC 1231 



3.1 The LHC experiments 

The four LHC experiments - i.e. the two large general-purpose 
and high-luminosity ATLAS and CMS detector systems as 
well as the heavy-ion-dedicated ALICE and the heavy-flavour- 
oriented LHCb experiments - have all detection capabilities in 
the forward direction very well adapted for the study of low-x 
QCD phenomena with hard processes in collisions with proton 
and ion beams: 

(i) Both CMS and ATLAS feature hadronic calorimeters in the 
range 3< |r|| <5 which allow them to measure jet cross- 
sections at very forward rapidities. Both experiments fea- 
ture also zero-degree calorimeters (ZDC, |r|| > 8.5 for neu- 
trals), which are a basic tool for neutron-tagging "ultra- 
peripheral" PbPb photoproduction interactions. CMS has 
an additional electromagnetic/hadronic calorimeter (CAS- 
TOR, 5.3< |r|| <6.7) and shares the interaction point with 
the TOTEM experiment providing two extra trackers at 
very forward rapidities (Tl, 3.1 < |r)| < 4.7, and T2, 5.5< 
|t|| <6.6) well-suited for DY measurements. 

(ii) The ALICE forward muon spectrometer at 2.5 < T| < 4, 
can reconstruct J/\\t and Y (as well as Z) in the dimuon 
channel, as well as statistically measure single inclusive 
heavy-quark production via semileptonic (muon) decays. 
ALICE counts also on ZDCs in both sides of the interac- 
tion point (IP) for forward neutron triggering of PbPb pho- 
toproduction processes. 

(iii) LHCb is a single-arm spectrometer covering rapidities 
1.8 < T| < 4.9, with very good particle identification capa- 
bilities designed to accurately reconstruct b and c mesons. 
The detector is also well-suited to measure jets, QQ and 
Z^ /j/j production in the forward hemisphere. 



3.2 Low-x QCD measurements at LHC 

Measurement at the LHC forward rapidities of any of the pro- 
cesses shown in Fig.|4]provides an excellent means to look for 
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signatures of high gluon density phenomena at low x. Four rep- 
resentative measurements are discussed in the last section of 
the paper. 



on the gluon structure and evolution at low-x. In the case of 
ALICE, the heavy-Q pp studies will have a natural extension 
in AA and pA collisions |41 1, providing a precise probe of non- 
linear effects in the nuclear wave-function. 



• Case study I: Forward (di)jets {pp, pA,AA) 

The jet measurements in pp collisions at Tevatron energies 
have provide valuable information on the proton PDFs (see 
Fig. |5] top). According to Eq. the measurement of jets with 
p T « 20 - 200 GeV/c in pp collisions at 14 TeV in the AT- 
LAS or CMS forward calorimeters (3< |r|| <5) allows one to 
probe the PDFs at x values as low as X2 ~ 10~ 4 — 10~ 5 . Fig. II II 
shows the actual log(xi,2) distribution of two partons scatter- 
ing at LHC and producing at least one forward jet as obtained 
with PYTHIA 6.403 |38|. As expected in forward scattering, 
the collision is very asymmetric with X2 (xi) peaked at ~ 10 -4 
(~ 1 0~ 1 ) and thus provides direct information on the low-x par- 
ton densities. Not only the single inclusive cross-section but the 
forward-backward dijet production, "Miiller-Navelet jets", is a 
particularly sensitive measure of BFKL 1391 as well as non- 
linear 1 6 1 parton evolutions. In the presence of low-x saturation 
effects, the Miiller-Navelet cross section for two jets separated 
by Ar| ~ 9 (and, thus, measurable in each one of the forward 
calorimeters) is expected to be suppressed by a factor of ^2 
compared to BFKL predictions |6|. A study is underway to de- 
termine the feasibility of both forward (di)jet measurements in 
CMS |30l . 
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PYTHIA 6.4: p+p-> je^+jet.,, xjs = 14 TeV 
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Fig. 11. log(;q,2) distribution of two partons colliding in pp collisions 
at v 7 ! = 14 TeV and producing at least one jet within ATLAS/CMS 
forward calorimeters acceptances as determined with PYTHIA 




1 year, pp, Vs = 14 TeV @ nominal lumin. 



— 100 1= 



a 

h- 



10 



ATLAS/CMS 
b 




-2 2 4 6 

r| of Q-hadron 

Fig. 12. Acceptances in (r\,pr) for open charm and bottom hadrons 
in the four LHC experiments for 1-year nominal luminosity 1151 



• Case study III: QQ photoproduction (electromagnetic AA) 

High-energy diffractive photoproduction of heavy vector 
mesons (J/\\i,Y) proceeds through colourless two-gluon ex- 
change (which subsequently couples to y — > QQ) and is thus 
a sensitive probe of the gluon densities at small x, see Eq. (@}. 
Ultra-peripheral interactions (UPCs) of high-energy heavy ions 
generate strong electromagnetic fields which help constrain 
the low-x behaviour of xG via quarkonia |42|, or other hard 
probes |43|, produced in y-nucleus collisions. Lead beams at 
2.75 TeV have Lorentz factors y = 2930 leading to maximum 
(equivalent) photon energies co max w y/R ~ 100 GeV, and cm. 
energies W™ ax « 160 GeV and W™ x w 1 TeV. From Eq. (0, 
the x values probed in yA — > J/\\i A processes at y = 2 can be 
as low as x ~ 10~ 5 . The ALICE, ATLAS and CMS experi- 
ment can measure 7/\|/,T — > e + e~ ,p + p~ produced in electro- 
magnetic PbPb collisions tagged with neutrons detected in the 
ZDCs (as done at RHIC |42|). Figure 13 shows the expected 
dimuon invariant mass distributions around the T mass pre- 
dicted by STARLIGHT 031 within the CMS central acceptance 
(T| <2.5) for an integrated PbPb luminosity of 0.5 nb~' |45 1. 



• Case study IV: Forward Drell-Yan pairs {pp, pA, AA) 



• Case study II: Forward heavy-quarks {pp, pA, AA) 

Studies of small-* effects on heavy flavour production at the 
LHC in two different approaches, based on collinear and kj 
factorization, including non-linear terms in the parton evolu- 
tion, lead to two different predictions (enhancement vs. sup- 
pression) for the measured c and b cross-sections |15|. The 
possibility of ALICE and LHCb to reconstruct D and B mesons 
in a large rapidity range (Fig. ll2l> will put stringent constraints 



High-mass Drell-Yan pair production at the very forward ra- 
pidities covered by LHCb and by the CMS CASTOR and 
TOTEM T2 detectors can probe the parton densities down to 
x ~ 10~ 6 at much higher virtualities M 2 than those accessible 
in other measurements discussed here. A study is currently un- 
derway in CMS 1 40 1 to combine the CASTOR electromagnetic 
energy measurement together with the good position resolution 
of T2 for charged tracks, to trigger on and reconstruct the e + e~ 
invariant mass in pp collisions at 14 TeV, and scrutinize xg in 
the M 2 and x plane. 
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Fig. 13. Expected invariant mass from yPb —* YPb* — > 

jU+^u" and yy — > /J + ,u~ processes predicted by Starlight 1 44 1 for 
UPC PWZ? collisions at v / 57m = 5.5 TeV in the CMS acceptance 



4 Conclusion 

We have reviewed the physics of non-linear QCD and high 
gluon densities at small fractional momenta x, with emphasis 
on the existing data at HERA (proton) and RHIC (nucleus) 
which support the existence of a parton saturation regime as 
described e.g. in the framework of the Colour Glass Conden- 
sate effective field theory. The future perspectives at the LHC 
have been presented, including the promising capabilities of the 
forward detectors of the ALICE, ATLAS, CMS and LHCb ex- 
periments to study the parton densities down to x ~ 10~ 6 with 
various hard probes (jets, quarkonia, heavy-quarks, Drell-Yan). 
The programme of investigating the dynamics of low-x QCD is 
not only appealing in its own right but it is an essential pre- 
requisite for predicting a large variety of hadron-, photon- and 
neutrino- scattering cross sections at very high energies. 
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